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Abstract 

Using Monte Carlo methods, the short-time dynamic scaling behaviour 
of two-dimensional critical XY systems is investigated. Our results for the 
XY model show that there exists universal scaling behaviour already in 
the short-time regime, but the values of the dynamic exponent z differ for 
different initial conditions. For the fully frustrated XY model, power law 
scaling behaviour is also observed in the short-time regime. However, a vi- 
olation of the standard scaling relation between the exponents is detected. 
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Recently much progress has been made in critical dynamics. It was discovered 
that universal scaling behaviour may emerge already in the macroscopic short- 
time regime [[1], [2|, [3], [|, |5|, [|. Extensive Monte Carlo simulations show that the 
short-time dynamic scaling is not only conceptually interesting but also practi- 
cally important, e.g. it leads to new ways for the determination of the critical 
exponents and the critical temperature || [7|, |8|, ||, [H| [11], |3| . 



For critical systems with second order phase transitions, comprehensive un- 
derstanding has been achieved. For a relaxational dynamic process of model A 
starting from an ordered state, the scaling form is given, e.g. for the k-th moment 
of the magnetization at the critical temperature, by || [| 

M^(t, L) = -Hd-2+v)/2 M (k)( b -z t ^ b -i L ^ 

where t is the time variable, L is the lattice size, t] and z represent the standard 
static and dynamic critical exponents. This scaling form looks similar to that 
in the long-time regime but it is now assumed to hold also in the macroscopic 
short-time regime after a microscopic time scale t m i C . 

For a relaxation process starting from a disordered state with small or zero 
initial magnetization, the scaling form for the fc-th moment at the critical tem- 
perature is found to be 

M (fc) (t,L,m ) = b' k{d - 2+ ^ l2 M {k \b- z t,b- l L,b X0 mQ). (2) 

Here it is important that a new independent critical exponent xq has been intro- 
duced to describe the dependence of the scaling behaviour on the initial magne- 
tization. The exponent Xq is the scaling dimension of the global magnetization 
M(t) and also of the magnetization density. Therefore, even if mo = 0, xq still 
enters observables related to the initial conditions. For example, for sufficiently 
large lattice size the auto-correlation has a power law behaviour 

A(t)^^(J2sms t (t)}-t- x (3) 

i 

with Si being a spin variable. The exponent A is related to xq by 
\ d a a ( d-2 + r]\ 1 

In general, for arbitrary initial magnetization mo between and 1 a generalized 
scaling form can be written down. 

Renormalization group calculations for the O(N) vector model show that the 
static exponents 7] and the dynamic exponent z in the scaling forms (0) and 
(|3[) take the same values as in equilibrium or in the long-time regime of the 
dynamic evolution where they are defined. This is also confirmed by Monte 
Carlo simulations for various critical magnetic systems with second order phase 
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transitions. Furthermore, numerical results indicate in good accuracy that the 
values of the exponents r\ and z are independent of the initial conditions, i.e. the 
values of the exponents in scaling form (0) are the same as in scaling forms (^j) and 
@. This prominent property has been used to extract all the dynamic and static 
exponents from the short-time dynamic scaling forms. Since the measurements 
are carried out in the short-time regime, the short-time dynamic approach is free 
of critical slowing down. One may believe that the scaling forms ([[]), (§) and 
(|3p hold also in the crossover regime where the time t is not small but also not 
asymptotically large. 

In this letter we investigate numerically whether critical systems with a Koster- 
litz-Thouless phase transition show also a clean scenario as stated above. For this 
purpose, we should answer the following questions: at the first, whether there 
exists universal scaling behaviour in the short-time regime; secondly, whether the 
dynamic exponent z and the static exponent r\ in the short-time scaling forms (|T|), 
(0) and (Q) are independent of the initial conditions and take the same values as 
in equilibrium or in the long-time regime where they are defined^; finally whether 
the scaling relations between the exponents for the second order phase transitions, 
e.g. that between xq and A in Eq. (^), hold also for Kosterlitz-Thouless phase 
transitions. 

As examples, let us consider the two-dimensional classical XY model and 
the fully frustrated XY (FFXY) model. The dynamic process starting from an 
ordered state has been investigated with Monte Carlo methods at the Kosterlitz- 
Thouless phase transition Tkt and below JT3] . Nice power law scaling behaviour 
was observed, and the static exponent r] extracted from scaling form (P is con- 
sistent with that measured in equilibrium. The dynamic exponent z is a constant 
for all temperatures at and below Tkt an d slightly smaller than 2.0, but within 



errors it coincides with the theoretical prediction z = 2 |14 . 

For the dynamic process starting from a disordered state, however, the situ- 
ation is complicated. For a quench to zero temperature, it was suggested some 
years ago that the dynamic scaling is violated [jT5fl . However, recent numerical 
simulations show that at least for a quench to not so low temperatures the mag- 
netization undergoes a universal power law initial increase if the initial magneti- 
zation is small but non-zero 0, [Hj]. Therefore, a comprehensive understanding 
of the short-time behaviour of this dynamic process is important and necessary. 

Our strategy is to simulate a dynamic relaxational process starting from a 
disordered state with zero magnetization and to measure the auto-correlation and 
the second moment of the magnetization. Then from the scaling forms (|2]) and (Q), 
we obtain the static exponent r\ and the dynamic exponent z and compare them 
with those measured from the scaling form ([!]). Indeed, to simulate a dynamic 

1 In principle, even if there exist the short-time scaling form as in Eqs. (|l]), (||) and (||), the 
dynamic exponent z and even the static exponents could be initial condition dependent and 
different from those in equilibrium since they are defined in different fixed points of the time t 
and different initial conditions. 
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process starting from a disordered state with zero magnetization is more difficult 
than that from an ordered state. This is very probably due to the effect of the 
vortices. We need very large lattices and update the system to fairly long time t. 

The XY model and the FFXY model in two dimensions can be defined by the 
Hamiltonian 

H = K J2 fa & ■ Sj , (5) 

<ij> 

where Si = (Si )X , Si tV ) is a planar unit vector at site i, and the sum extends over 
the nearest neighbours. In our notation, K is just the inverse temperature. Here 
fij take the values +1 or —1, depending on the model. For the XY model, /y = 1 
on all links. A simple realization of the FFXY model is by taking = — 1 on 
half of the vertical links (negative links) and +1 on the others (positive links) 



It 



In this paper, only the dynamics of model A is concerned. The dynamics of 
model A is a relaxational dynamics without energy and magnetization conserva- 
tion. Starting from a completely random configuration, we update the system at 
the temperature Tkt or below with the standard Metropolis algorithm. The trial 
state of each spin is randomly taken in the unit circle since the acceptance rate 
is high in the short-time regime of the dynamic evolution. Updating is stopped 
at Monte Carlo time step t = 1500. The procedure is repeated with another 
initial configuration and different random numbers. We have tested that for up- 
dating time t = 1500, a big lattice size like L = 512 is needed. An updating 
time t = 1500 is necessary for reliable measurements of the observables since the 
microscopic time scale is here around t m j C ~ 100 or more. For the average a total 
of 800 samples has been taken. Errors are estimated by dividing the samples into 
four groups. 

We measure the auto-correlation 

Mt) = ji(E3 i {o)-§ i {t)) (6) 

i 

and the second moment of the magnetization 



(7) 



For the FFXY model, the second moment must be calculated in four sublattices 
separately, and the results are summed up after average. Keeping in mind that 
spatial correlation is very small in the short-time regime of the dynamic evolution, 
from the finite size scaling in Eq. (0) one may easily deduce the short-time power 
law behaviour for the second moment 

M {2 \t)~t y (8) 



3 



with 

y = (2- V )/z. (9) 

In Fig. [I], the auto-correlation for both the XY and the FFXY model is dis- 
played in log- log scale for different temperatures at and below T KT . After a cer- 
tain microscopic time scale t m i C , power law behaviour is clearly observed. Careful 
analysis shows that for temperature T not so far from Tkt, tmic ~ 100 — 200 but 
as T decreases, t mic gradually increases. Actually, this tendency can also be seen 
roughly by eyes from the figure. For example, for the XY model at T = 0.70, 
tmic ~ 300 — 400. This increase of t mic for lower temperatures was also noticed 
in the measurements of the critical initial increase of the magnetization [pi], |i~7H . 



This phenomenon is somehow understandable since at low temperatures the con- 
figuration tends to be frozen. As the temperature T approaches zero, we may 
believe that t mic becomes considerable big. Numerical simulations of the short- 
time behaviour for very low temperatures are much more difficult. It requests 
long updating time t and correspondingly large lattices. This might explain the 
violation of dynamic scaling observed for a quench to zero temperature Jl5| 



From the slopes of the curves in Fig. |T|, we measure the exponent A. The 
results are given in Tables [I] and |2|. As the temperature decreases, A becomes 
smaller. We should mention that the errors here are only those estimated by 
dividing the total samples into four groups. Other errors as the fluctuation in 
time direction and a possible remaining effect of t m i C have not been taken into 
account^. In equilibrium, it is well known that critical slowing down is much 
more severe in the XY and the FFXY model than in the Ising model. Numerical 
measurements in the XY systems are very difficult. However, the situation is 
completely different for short-time dynamic measurements. For a lattice size 
L = 512 and updating time t = 1500, we already obtain rather accurate results 
with a total number of 800 samples. 

In Fig. H, the second moment for both the XY and FFXY model is displayed 
in log- log scale for different temperatures at and below Tkt- Again we observe 
power law behaviour even though the fluctuations here are apparently larger 
than those of the auto-correlation. From the slopes of the curves, we measure the 
exponent y. The results are also listed in Tables [l] and |[ Interestingly, for the 
FFXY model the exponent y obviously first increases and then decreases as the 
temperature goes down. This might be related to the rapid drop of the exponent 
A near Tkt- 

In Tables [I] and |2], the static exponent r\ is estimated from a dynamic process 
starting from an ordered state fT3|]. The results are in good agreement with and 
improve the measurements in equilibrium. For the XY model, taking the expo- 
nents 6 and r\ obtained in Refs. ]16], |T3j| as input, one can compute separately 
the dynamic exponent z from A and y for different temperatures. The results are 

2 The effect of t m ic may sometimes be a kind of correction to the scaling which may not 
disappear suddenly, e.g. when the correction obeys power law. 
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given in Table |I[ The values of z obtained from the auto-correlation and the sec- 
ond moment agree well within statistical errors. This fact strongly supports the 
scaling forms in Eqs. (fj) and (|3|) and the scaling relations between the exponents 
in Eq. (|J) and (||]). However, the dynamic exponent z is clearly different from that 
measured from a dynamic process quenched from an ordered initial state, which 
is denoted by Z\ in Table |l|. As the temperature decreases, Z\ remains a constant 
near 2, but z is apparently bigger than 2 and increases gradually. The differ- 
ence is 15 to 20 percent and we regard it as prominent. This indicates that the 
dynamic scaling behaviour of the XY systems with a Kosterlitz-Thouless phase 
transition is more complicated than the simple scenario of the dynamic scaling 
for critical systems with second order phase transitions. The dynamic XY model 
has a non-trivial crossover behaviour for the time t macroscopically not small but 
also not asymptotically large. As mentioned before, the microscopic time scale 
t mic grows as the temperature decreases. When the temperature approaches zero, 
it might happen that the system evolves already into the crossover regime at the 
time t ~ t mic , and thus the simple scaling forms (Q) and (|3|) are not observed. 

For the FFXY model, the situation is even more complicated. Taking r] as 
input, the resulting dynamic exponent z from the second moment interestingly 
shows a non-monotonous dependence on the temperature. This was also qualita- 



tively observed by other authors in calculating the domain growth |19| . However, 
if we take the exponent 9 as input, the dynamic exponent estimated from the 
auto-correlation with the scaling relation (J|) does not coincide with that from 
the second moment. The difference can not be explained by the statistical er- 
rors. Since the exponent y from the second moment is not directly related to 
the initial conditions, we believe that the scaling relation (|]) holds in any case, 
i.e. the exponent z computed from y should be correct. However, for the FFXY 
model the scaling relation (§) is violated. The violation can be described by 
5 = A — dj z + 9. Values for this quantity are given in Table 0. The origin of 
this violation may probably be traced back to the chiral degree of freedom. Since 
the chiral transition temperature T c is above T^t-, ordering dynamics of the chi- 
ral magnetization may affect the behaviour of the auto-correlation. However, a 
complete understanding of this problem remains open. 

In conclusions, we have numerically investigated the short-time dynamic be- 
haviour of the two-dimensional XY and the fully frustrated XY model at the 
temperature T KT and below. Our results show that there exists indeed dynamic 
scaling in the short-time regime. However, the dynamic exponent z may depend 
on the initial conditions. It indicates that a non-trivial crossover behaviour occurs 
when the dynamic evolution crosses over from macroscopic short-time regime to 
the long-time regime. This scenario is very different from that of critical systems 
with second order phase transitions. 

Acknowledgements: Work supported in part by the Deutsche Forschungs- 
gemeinschaft; Schu 95/9-1 and SFB 418. 
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T 0.90 0.86 0.80 0.70 



A 




0.625(4) 


0.600(3) 


0.569(2) 


0.552(4) 


y 




0.766(16) 


0.775(15) 


0.780(19) 


0.773(24) 


e 




0.250(1) 


0.264(5) 


0.290(1) 


0.287(3) 


V 




0.244(5) 


0.212(4) 


0.178(2) 


0.143(3) 


z = (2- 


v)/y 


2.292(48) 


2.307(45) 


2.335(57) 


2.402(75) 


z = d/{9 


+ A) 


2.286(11) 


2.315(16) 


2.328(06) 


2.384(14) 






1.96(4) 


1.98(4) 


1.94(2) 


1.98(4) 



Table 1: The exponents A and y measured for the XY model. Values of the 



exponent 8 are taken from Ref. fl6 |, while rj and z% are from Ref. fl3 |. z\ is the 



dynamic exponent z measured from a quench from an ordered initial state [13 



T 


0.44 


0.40 


0.35 


0.30 


A 


0.808(3) 


0.622(2) 


0.577(2) 


0.561(2) 


y 


0.747(14) 


0.835(20) 


0.809(15) 


0.801(13) 


6 


0.079(4) 


0.181(5) 


0.245(3) 


0.263(2) 


V 


0.243(4) 


0.140(2) 


0.107(1) 


0.086(2) 


z = ( 2 ~v)/y 


2.352(44) 


2.228(53) 


2.340(43) 


2.390(39) 


s = x-d/z + e 


0.037(17) 


-0.095(22) 


-0.033(16) 


-0.013(14) 


Zl 


1.93(4) 


1.95(3) 


1.99(7) 


1.97(3) 



Table 2: The exponents A and y measured for the FFXY model. Values of the 



exponents 9 are taken from Ref. fi~?|| , while rj and z\ are from Ref. |I7]]. z\ is the 



dynamic exponent z measured from a quench from an ordered initial state [13|. 
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Figure 1: Auto-correlation in log-log scale (a) for the XY model at the tempera- 
tures T = 0.90, 0.86, 0.80 and 0.70 (from below); (b) for the FFXY model at the 
temperatures T = 0.44, 0.40, 0.35 and 0.30 (from below). 
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Figure 2: Second moment in log-log scale (a) for the XY model with the solid, 
dotted, dashed and long dashed lines for the temperatures T = 0.90, 0.86, 0.80 
and 0.70 respectively; (b) for the FFXY model with the solid, dotted, dashed and 
long dashed lines for the temperatures T = 0.44, 0.40, 0.35 and 0.30 respectively. 
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